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ABSTRACT 

An  analytical  model  of  a  jet  injected  normally  from  a  flat 
plate  into  a  uniform  crossing  flow  was  formulated  to  provide 
a  simplified  method  of  predicting  the  interference  effects 
arising  from  the  complex  flow  fields  induced  by  ship  bow 
thrusters.   This  model  was  an  extension  of  previous  work 
based  upon  a  description  of  the  jet  as  a  series  of  distributed 
vortices.   The  analysis  takes  into  account  the  position  of 
the  effective  source  of  the  jet  and  the  blockage  due  to  the 
presence  of  the  jet  in  the  crossflow.   For  representative 
jet-to-crossf low  velocity  ratios,  the  flow  field  and  pressure 
distributions  were  calculated  utilizing  different  combinations 
of  effective  source  position  and  blockage.   The  accuracy  of 
the  model  was  evaluated  by  comparison  with  the  available 
experimental  data.   Although  good  agreement  was  achieved  for 
large  portions  of  the  interaction  field,  several  regions  were 
identified  requiring  further  analytical  description. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION  9 

II.  BACKGROUND 13 

III.  REVIEW  OF  WOOLER'S  VORTEX  MODEL  18 

IV.  EXTENSION  OF  WOOLER  MODEL  24 

V.  IMPROVEMENT  OF  CONFORMING  VORTEX  MODEL  29 

VI.  CURRENT  MODEL  CAPABILITY  32 

VII.  CONCLUSIONS 33 

VIII.  RECOMMENDATIONS 34 

APPENDIX  A   FORMULATION  OF  INTERFERENCE  VELOCITY 

COMPONENTS 51 

COMPUTER  PROGRAM  CONFORMING  VORTEX  MODEL  62 

LIST  OF  REFERENCES 68 

INITIAL  DISTRIBUTION  LIST  70 


LIST  CF  FIGURES 

1.  Sketch  Of  Bow  Thruster-Crossf low  Interaction 35 

2.  Description  Of  Coordinate  Systems 36 

3.  Wooler  Horseshoe  Vortex  System 37 

4.  Comparison  Of  Pressure  Coefficient  Contours 

Velocity  Ratio  =  4 38 

5.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  8 39 

y 

6.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  11.3 40 

7.  Conforming  Vortex  Loop  System  41 

8.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  4.0 42 

9.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  8.0 43 

10.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  ,11.3 44 

11.  Comparison  Of  Wooler  (a)  And  Conforming 
Vortex  taodel  (b)  Effective  Radii  For 

Induced  Velocities 45 

12.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =8.0 46 

y 

13.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =4.0 47 

14.  Comparison  Of  Pressure  Coefficient  Contours 

U./U  =  8.0 48 

15.  Comparison  Of  Pressure  Coefficient  Contours 

Uj/U  =  11.3 49 

16.  Variation  In  The  Induced  Force  Due  To  Jet-Cross flow 
Interaction — — — 50 


Al.    Coordinate  Systems  For  Conforming  Vortex  Model  59 

A2.    Geometry  Of  Bound  Vortex  Integration  60 

A3.    Geometry  Of  Trailing  Vortex  Integration  61 


ACKNOWLEDGEMENT 

The  author  wishes  to  express  his  appreciation  to  Associate 
Professor  Robert  Nunn  for  his  advice  and  guidance  during  the 
development  of  this  model.   Without  his  valuable  assistance, 
this  could  not  have  been  completed. 

Many  thanks  to  my  wife,  Helga,  for  her  infinite  under- 
standing and  patience  during  the  long  hours  needed  for  its 
completion. 


I.  INTRODUCTION 

In  recent  years  there  has  been  an  increase  in  the  instal- 
lation and  use  of  hull  mounted,  tunneled  bow  thrusters  to 
improve  ship  maneuverability  in  restricted  situations,  such 
as  dynamic  position  holding  and  confined  waters  docking. 
This  increase  has  created  renewed  interest  in  the  development 
of  an  ability  to  predict  bow  thruster  performance.   To  date, 
there  has  been  relatively  little  pre-installation  design 
consideration  regarding  precise  power  requirements  or  ultimate 
operational  performance.   The  limited  testing  that  has  been 
conducted  /~l-47*  has  resulted  in  a  recurring  major  problem, 
namely,  radical  variations  in  effective  turning  moment  with 
minor  changes  in  ship  speed.   Chislett  and  Bjorheden  /— 4  7 
noted  that  "an  area  of  low  pressure  is  created  downstream 
of  the  discharging  jet  producing  a  resultant  suction  force 
that  has  a  shifting  center  of  action  with  changes  in  ship 
speed,"  thereby  altering  the  effective  moment  produced  by 
the  thruster.   The  ability  to  predict  the  magnitude  of  this 
resultant  suction  force  and  its  center  of  action  would  permit 
designers  to  accurately  determine  powering  and  control  re- 
quirements for  thrusters. 

This  problem  of  radical  variations  in  effective  side 
forces  and  turning  moments  with  only  minor  changes  in  ship 
speed  may  be  better  understood  by  considering  the  general 
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equations  of  motion  of  a  ship  moving  in  the  horizontal 
plane  of  the  sea.   From  /""5_7  these  equations  are: 

m(u  -  rv)  =  X 

m(v  -  ru)  =  Y 

I  r  =  N 

z 

where  m  =  mass  of  the  ship 

u  =  velocity  in  the  x-direction 

v  =  velocity  in  the  y-direction 

r  =  angular  velocity 

I  =  mass  moment  of  inertia 

X  =  excitation  force (s)  in  x-direction 

Y  -  excitation  force (s)  in  y-direction 

N  =  excitation  torque  (s) 

(dots  indicate  differentiation  with  respect  to  time) 
To  minimize  the  complexity  of  the  following  brief  analysis, 
all  forces  and  torques  will  be  referred  to  the  shin's  center 
of  gravity.   See  Fig.  1. 

When  these  equations  are  transformed  into  their  linearized 
version  /~~5_7r  they  become 

-  X  (u  -  u  )  +  (m  -  X»)  u   =  X 

u  u   uo  u   o 

-YV+  (m-Y')v-  (Y   -mu)r-Y*r=Y 
v  v       r     o      r 

-  N  v  -  N»v  -  N  r  +  (I   -  N«)r  =  N 

v     v     r      z    r' 

... 

where  the  subscripts  u,  u,  v,  v,  r,  r,  denote  differentiation 

with  respect  to  the  variable  subscripts,  and  u   is  constant 
reference  surge  velocity.   The  left  sides  of  these  equations 
explicitly  represent  the  velocity  and  acceleration  dependent 
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reactions  due  to  drag,  added  mass,  added  inertia  and  other 

coupled  motions;  while  the  right  sides  implicitly  represent 

excitation  force (s)  and  moment (s) ,  which  are  in  this  analysis 

due  to  the  thruster.   Since  this  is  a  linearized  formulation 

utilizing  small  angle  perturbations,  the  excitation  force  in 

the  x-direction  is  negligible  compared  to  that  in  the 

y-direction: 

X  <<  Y 

therefore  X=0.   The  y-direction  excitation  force  can  be 

expressed  as 

Y  =  Y.  +  Y.  (1) 

3     1 

where  Y.  =  m.U.,  the  ideal  jet  thrust  produced  by  the 
thruster,  and  Y.  is  the  induced  force  created  by  the  inter- 
action of  the  thruster  discharge  with  the  flow  past  the  ship. 
Utilizing  experimental  data  from  /~"4  7  in  the  range, 

0.1  £  ~  <  0.25,  Y.  can  be  approximated  by  a  linear  function 

J 
of  ship  speed 

Y. 

~  =  0.22  +  k(H_)  (2) 

J  J 

where  -3.0  <  k  <  -1.5.   Substituting  a  nominal  value  of 
k  =  -2.0  and  Y.  =  m.U.  into  equation  (2)  gives 


Y.  =  m.U.[o.22  -  2.0  (£-)  ] 


(3) 


Substituting  equation  (3)  into  equation  (1)  results  in 
Y  =  i^U..  [l.22  -  2.0(~)  ] 
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which  explicitly  couples  the  y-direction  motion  to  the 
x-direction  motion,  normally  uncoupled  in  the  linear 
approximation.   In  addition  to  the  coupling,  in  the  relevant 

velocity  ratio  range,  0.1  <  rr—  <    0.25,  there  is  an  appreciable 

J 

reduction  in  the  net  thruster  force  due  to  motion  in  the 

x-direction.   Furthermore,  the  moment  produced  by  the  thruster 
is  a  function  of  the  net  effective  thrust,  Y,  and  a  character- 


istic moment  arm,  x .  : 


■   =  XtY 


As  reported  by  Refs.  /"l-3/,  the  degradation  of  net  effective 
thrust  is  accompanied  by  a  shift  in  the  center  of  action  of 
the  induced  force,  Y.,  which  results  in  a  change  in  the  char- 
acteristic moment  arm.   This  shift  in  the  thruster' s  net 
effective  thrust  center  of  action  is  a  direct  consequence  of 
the  coupling  of  x-direction  motion  with  turning  moment. 

This  Simplified  analysis,  which  neglected  all  excitations 

other  than  the  thruster,  points  out  the  complex,  coupled  mo- 
I 

tions  arising  from  the  interaction  of  the  thruster  discharge 
and  the  flow  past  the  ship.   To  fully  understand  this  highly 
non-linear,  three-dimensional  flow  situation  extensive 
analytical  and  experimental  work  is  required.   This  thesis 
will  attempt  to  provide  some  understanding  of  the  interaction 
of  a  jet  discharging  into  a  uniform  cross-flow  in  order  to 
develop  some  ability  to  predict  the  degradation  of  net  effective 
thrust. 
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II.  BACKGROUND 

A  survey  of  the  literature  pertaining  to  the  installation 
and  operation  of  bow  thrusters  revealed  a  surprising  lack  of 
experimental  data,  considering  the  number  of  units  currently 
installed.   English  /T,2_7 \    Stuntz  and  Taylor  /~"3_7,  and 
Chislett  and  Bjorheden  /"~4_7  present  some  data  from  model  tests 
conducted  at  various  facilities,  but  in  these  works  there  is 
no  formalized  recording  of  full  scale  test  or  actual  instal- 
lation test  data.   However,  further  investigation  led  to  a 
related  field:  Vertical  or  Short  Take-Off  and  Landing  Aircraft 
(V/STOL)  research.   The  concept  of  using  air  jets  positioned 
in  aircraft  wings  to  provide  vertical  take-off  and  landing 
capability  correlates  extremely  well  with  flow  patterns  pro- 
duced by  a  bow  thruster.   As  reported  by  Bradbury  and  Wood 
J"~6  7  for  incompressible  flow,  the  jet  (thruster  discharge) 

path  and  induced  flow  are  dependent  mainly  on  the  momentum 

2    2 

flux  ratio,  p.U.  /pU  ,  and  independent  of  Reynolds  number. 

Margason  /"7  7  and  Gordier  /"~8  7  also  reported  that  the 
effective  velocity  ratio  is  the  predominate  characteristic 
in  determining  the  path  of  the  jet.   Therefore,  until  more 
extensive  ship  or  model  data  have  been  collected,  the  re- 
sults of  V/STOL  research  are  used  as  a  basis  for  predictive 
theories. 

Although  there  has  been  a  good  deal  of  research  conducted 
in  the  area  of  air  jets  in  cross-flows,  the  majority  of  the 
work  has  been  concerned  with  analytically  or  empirically 
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defining  the  jet  trajectory  and  jet  cross-section  geometry. 
The  interference  effects  of  the  interacting  flows  have  re- 
ceived relatively  little  attention.   A  fluid  jet  injected 
into  a  crossing  stream  has  previously  been  described  as  a 
turbulent,  three-dimensional,  highly  non-linear  flow  problem 
which,  even  with  considerable  simplification,  requires  ex- 
tensive computer  time  to  numerically  solve  the  appropriate 
form  of  the  Navier-Stokes  equations  /~"9_7.   Therefore,  early 
investigators  have  sought  simplified  methods  for  determining 
the  interaction  effects,  such  as  pressure  and  velocity  dis- 
tributions, and  have  realized  the  necessity  of  knowing  the 
geometry  of  the  jet:   trajectory  and  shape. 

Experimental  data  form  the  basis  for  most  trajectory 
formulations.   Either  pressure  or  velocity  measurements  have 
been  taken  in  the  flow  field  encompassing  the  jet-cross flow 
interaction  region  and  curves  fit  through  points  of  maximum 
pressure  or  velocity.   These  curves  were  then  compared  with 
flow  visualizations.   In  some  instances,  semi-empirical  for- 
mulations were  derived  from  the  conservation  laws  combined 
with  experimentally  determined  constants.   Abramovich  /~~107 
qualitatively  described  the  turbulent  jet  in  a  deflecting 
flow  and  presented  some  empirical  methods  for  predicting 
the  trajectory  of  the  deflected  jet.   Jordinson  /Tl  /  pre- 
sented trajectory  data  by  recording  and  plotting  contours 
of  total  pressure  coefficients.   Keffer  and  Baines  /"~12_7 
presented  experimental  results  from  which  it  was  determined 
that  for  various  jet  strengths,  the  jet  trajectories  could 
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be  represented  by  a  single  function.   Gordier  ,  —  3_7,  the 
only  reported  source  of  water  jet-water  crossflow  experiments, 
presented  an  empirical  jet  trajectory  formulation  based  on 
a  curve  passed  through  experimental  points  of  maximum  total 
pressure.   Sucec  and  Bowley  /~~9_7  formulated  an  analytical 
expression  for  the  jet  trajectory  utilizing  previous  experi- 
mental information  and  the  assumption  that  the  distributed 
pressure  force  and  entrained  momentum  flux  could  be  approxi- 
mated by  an  aerodynamic  drag  force.   Margason  /"" ij   utilized 
flow  visualization  techniques  and  pressure  measurements  to 
formulate  a  trajectory  equation.   In  addition,  his  paper  also 
reviewed  and  compared  the  results  of  a  number  of  other 
studies  concerned  with  trajectory  prediction.   All  the  above 
predictive  formulations  produced  trajectories  within  a  range 
of  uncertainty  that  can  be  attributed  to  expected  experimental 
error,  differences  in  test  procedures  and  facilities,  and  in 
the  case  of  analytical  formulations,  simplifying  assumptions. 
In  view  of  this,  it  was  concluded  that  any  one  of  them  re- 
flected the  current  state  of  the  art  in  trajectory  prediction. 

The  ability  to  predict  the  jet  trajectory  is  of  definite 
importance,  but  the  primary  objective  of  this  jet-cross flow 
research  is  to  provide  insight  into  the  changes,  in  the  other- 
wise uniform  flow,  created  by  injecting  a  jet.   Experiments 
conducted  by  Bradbury  and  Wood  /"~6_7»  Vogler  /~~13_7/  Fearn 
and  Weston  /""  14__7,  McMahon  and  Mosher  /""  15_7#  Kamotani  and 
Greber  /~~16  7  and  others,  provide  measured  values  of  the 
pressure  distribution  on  the  surface  surrounding  the  jet 
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orifice  which  have  been  used  to  compare  the  accuracy  of 
predictive  models.   During  the  past  few  years,  a  number  of 
approximate  methods  to  oredict  the  pressure  distribution  has 
been  developed:   Wooler  /T7 ,187 ,    Wu  and  Wright  /197/  Adler 
and  Baron  /~~207  and  Schmitt  /~~ 217.   Again,  because  of  the  com- 
plex nature  of  the  flow,  all  these  models  have  utilized  some 
degree  of  empiricism.   Either  the  model  is  formulated  using 
experimentally  obtained  trajectory  expressions  or  the  model 
is  formulated  using  analytical  expressions  whose  coefficients 
are  selected  to  provide  best  fit  with  experimental  data. 
The  majority  of  the  models  are  based  on  integral  techniques 
with  simplifying  assumptions,  such  as  (i)  the  representation 
of  the  entire  flow  as  two-dimensional  /197,  (ii)  the  external 
flow  is  irrotaticnal,  uncompressible  and  steady-state  /~20j   or 
(iii)  the  flow  is  inviscid  except  that  viscosity  is  the 
mechanism  that  leads  to  entrainment  /T8  7» 

With  the  aid  of  these  assumptions,  the  following  approaches 
were  made.   Adler  and  Baron  formulated  their  model  by  integrat- 
ing momentum  equations  to  describe  the  jet  mixing  field,  without 
using  empirical  trajectory  data.   However,  numerous  parameters 
were  derived  using  empirical  correlations  from  prior  research. 
Their  results  provide  satisfactory  agreement  with  experiment, 
but  from  the  quantity  of  empiricism  involved,  one  might  ex- 
pect such  agreement.   Wooler  /"~177  followed  by  Wu  and  Wright 
/  19__7  utilized  blockage-sink  representations  where  entrain- 
ment of  crossflow  fluid  was  handled  analytically  with  suitably 
chosen  coefficients.   Again  the  results  had  satisfactory 
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agreement  with  experiment.   However,  an  earlier  effort  by 
Wooler  /"" 18  7  minimized  empiricism  to  the  extent  that  only 
an  experimentally  determined  jet  trajectory  was  necessary 
to  complete  the  model.   By  thus  restricting  the  number  of 
empirical  parameters,  this  model  has  considerable  appeal 
from  an  engineering  point  of  view.   Since  visualization  of 
actual  jet-crossf low  situations  has  shown  that  the  jet  is 
deflected  in  the  crossflow  direction  and  forms  two  contra- 
rotating  trailing  vortices,  this  interference  model  was 
based  on  the  representation  of  the  jet  by  a  distribution  of 
vorticity.   Incorporating  arguments  from  the  aerodynamic 
theory  of  lifting  bodies  this  vorticity  distribution  was 
quantitatively  described  along  the  experimentally  determined 
trajectory.   By  virtue  of  minimal  empiricism  and  therefore 
greater  appeal,  this  model  was  selected  as  the  basis  for 
continued  investigation  in  developing  a  better  predictive 
tool. 
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III.  REVIEW  OF  WOOLER'S  VORTEX  MODEL 

Before  any  attempt  to  develop  an  improved  model  was 
made,  a  detailed  re-derivation  of  the  basis  formulation  was 
performed  in  order  to:   (1)  gain  a  better  understanding  of 
the  logic  of  formulation  and  (2)  confirm  the  reported  results 
A  review  of  the  distributed  vortex  method  follows. 

One  area  of  interest  was  the  geometry  of  the  problem. 
The  empirical  trajectory  equation  used  was 


=   B  [cosh(z/Bd)   -  l] 


x/d   =   B  |cosh(z/Bd)   -  lj  (4) 

where  x  is  measured  in  the  direction  of  the  mainstream  and 
z  is  measured  in  the  direction  of  the  exiting   jet.   The 

coefficient  B  was  empirically  determined  from  Jordinson  /Tl7 

2 

and  is  equal  to  0.19  (U./U)  ,  where  U.  is  the  jet  velocity 

and  U  is  the  mainstream  velocity.   To  further  describe  the 
problem  i  (see  Fig.  2),  a  system  of  natural  coordinates 
attached  to  the  jet  is  adopted.   Relative  to  the  fixed 
coordinate  system  and  an  arbitrary  point  jx  ,y  1  on  the 
surface  surrounding  the  jet  orifice,  the  natural  coordinates 
are  given  by 


E,   =      -    [zsina  +  (x~xp*  oosct  J 


n  -  Yp  (5) 

C  =   zcosa  -  (x-x  )  sina 

P 

where  £  is  the  coordinate  tangent  to  the  jet;   S,  the  coordi- 
nate normal  to  the  jet  in  the  direction  of  the  center  of 
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curvature;  r\ ,    the  coordinate  perpendicular  to  £  and  £  ; 
and  a  ,  the  angle  between  the  x  and  £,   directions. 

Another  extremely  important  part  of  Wooler's  formulation 
was  the  calculation  of  the  distribution  of  vorticity  within 
the  jet.   This  distribution  of  vorticity  was  determined  by 
adopting  methods  commonly  found  in  the  aerodynamic  theory  of 
lifting  surfaces  /"~227.   According  to  this  theory,  regions  of 
flow  external  to  the  jet  are  taken  to  be  irrotational  and  the 
deflection  of  the  jet  is  due  to  a  purely  inviscid  mechanism. 
This  inviscid  mechanism  is  expressed  as  a  balance  between  the 
forces  due  to  pressure  differences  across  an  element  of  the 
jet  and  the  centrifugal  forces  associated  with  the  jet  curva- 
ture.  These  pressure  differences  define  corresponding  velocity 
differences  according  to  Bernoulli's  Theorem.   These  resulting 
velocity  differences,  in  turn,  are  related  to  circulation 
according  to  Kelvin's  definition  of  circulation.   Therefore, 
with  the  above  assumptions,  the  circulation  is  related  to 
jet  curvature,  and  this  relationship  is 

2 


r  _  -  i  (  uj 


U     T   R\U 


(6) 


where  T   is  the  circulation  per  unit  length,  ds,  along  the 
jet;  R,  the  local  radius  of  curvature  of  the  jet;  and  d,  the 
diameter  of  the  jet  orifice.   Utilizing  the  empirical  tra- 
jectory expression  /~~Eq.  (427  in  the  arc  length  derivative 
results  in 


ds   _  1  dz  1-jk 

R      Bd  cosh(z/Bd) 
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Incorporating  Eq.  (7)  into  Eq.  (6)  and  integrating  over  an 
element  of  jet  extending  from  z,  to  z»  results  in  the  follow- 
ing expression  for  the  circulation  of  a  jet  element 


K      -      JudQ 


2   „ 


-1         Z2  -1         Zl 

tan        e        /Bd  -   tan        e        /Bd 


(8) 


In  order  to  determine  the  flow  field  disturbance  on  the  plate, 
the  jet  was  divided  into  a  number  of  vortex  elements  whose 
strengths  were  calculated  using  Eq.  (8).   These  elements  are 
the  origins  of  characteristic  horseshoe  vortex  systems  where 
the  cross  member  is  bound  in  the  jet  and  the  trailing  sides 
are,  according  to  Wooler's  assumption,  tangent  to  the  jet  and 
separated  by  a  distance,  d,  the  diameter  of  the  jet  orifice 
(Fig.  3) .   Errors  associated  with  taking  the  trailing  vortices 
tangent  to  the  jet  instead  of  containing  them  within  the  jet 
as  in  the  actual  flow  situation  were  assumed  by  Wooler  to  have 
negligible  effect  on  the  velocity  field  along  the  plate.   Each 
of  these  horseshoe  vortex  systems  produces  an  interference 
velocity  on  the  surface  surrounding  the  jet  orifice.   To  de- 
termine this  interference  velocity  at  any  arbitrary  point  on 
the  surface,  the  Law  of  Biot  and  Savart  /"~237 


q    =  s= 


K   /  ds  x  r 


L 
is  applied  to  each  bound  vortex  and  its  associated  trailing 

vortices.   For  the  bound  vortex,  -   /2  <  L  <  /2 ,  and  for  the 

trailing  vortices,  z^   <   L   <   <x>    ,   were  the  integration  limits. 

These  integrations  produced  closed-form  solutions  for  the 

interference  velocities  due  to  individual  jet  elements.   The 
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total  interference  velocity  components  due  to  the  entire  jet 
are  determined  as  the  sum  of  elemental  contributions  and 
are  given  by 


2  N   r 


8    *vu v  Z, 


tan 


v 

U 


3VU  i 


i=l 
2    N  ,_ 

i=l  - 


z . , , /Bd    2 . /Bd 
l  +  l        l 

-  i,  e -  e 

.  (z.^,  +  z.  )/Bdy 
.+e   l+l    i  ' 


x  (-w.+w^-w,) sina. +u,cosa.l 


(-v2+v3) tan 


z.+1/Bd    zi/Bd 

-  u   -3        -  e 


-1  +  e 


(zi+l  +  Zi)/Bd 


where  N  is  the  number  of  jet  elements  and  z.  and  z. ,,  are  the 

J  l      l+l 

endpoints  of  the  general  element.   The  parameters  u,  v  and 
w  are  defined  as  geometric  coefficients  resulting  from  the 
integration  along  the  jet.   The  subscripts  1,  2,  3  refer  to 
the  bound,  left  and  right  (as  viewed  from  upstream)  vortex 
contributions  respectively.   They  are  given  by 


u,  =  -£ 


1      e2- 


g  +  is 


-  H 


+C        lyt2+z2+    (n+J5)2        Y^2+c2+    (n-  h) 


v-,  = 


2      ^T,      ZT2 
c  +(n-  **) 


l  + 


Jz2+z2+   <n-  ^)2 


v.,   = 


3  "  "2 


C  +(n+  *) 


7       + 


U/s2+  ^2+  (^  +  ?5)2. 
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wl = 


FT 


?w 


n  +  h 


-  h 


1  +   z2  +    (n   +  H)2 


K2  +  t,2  +    (n  -  h)2 


w2  =  -n 


z,2  +  (  n-  is)2 


'V 


C2  +  C2  +    (n  -  h)2 


r\   +   H 
w3   =   


C2  +    (n  +  ^s)2 


'V 


S2  +  c2  +   (n  +  ^)2 


and  5f  C  ,  and  n  are  given  by  Eq.  (5)  for  each  point  on  the 
surface.   After  calculation  of  the  interference  velocity 
vector  at  a  sufficient  number  of  Doints  on  the  surface,  the 
pressure  coefficinet  at  each  point,  C  ,  is  given  by 


C  = 


P   h 


P  -  P 

00 

i     2,  2 

j-        =  1  -  q  /U 


whe 


%    pU 
re  q2/U2  =  fU  ±   2u J   +(-§^)  •   The  double  u  and  v 


components  are  the  result  of  the  use  of  an  image  system  to 
establish  the  surface  as  a  solid  boundary  thereby  negating 
the  w  component  of  velocity  in  this  plane.   Upon  substitution 
and  simplification,  one  has 


Cp  =  -4 


&♦(&)    +® 


(9) 


After  completion  of  the  review  and  the  incorporation 
of  several  necessary  corrections  to  the  published  work  /T8__7> 
the  entire  formulation  was  coded  for  computer  solution.   The 
resulting  plots  of  pressure  contours  for  representative 
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velocity  ratios  (Figs.  4,5,6)  show  fair  agreement  with  the 
experiments  of  Bradbury  and  Wood  ,  6  /.   However,  these 
plots  also  point  out  serious  deficiencies  in  the  model.  In 
the  far  field,  in  the  arc  35   <  9  <  145   (8  being  measured 
counter  clockwise  from  the  ray  extending  downstream  from 
the  jet  origin) ,  the  method  represents  the  jet  interference 
effects  well,  but  in  both  the  upstream  and  downstream 
portions  of  the  surface,  there  is  a  definite  lack  of  agree- 
ment.  In  the  upstream  area,  the  blockage  effect  due  to  the 
presence  of  the  jet  in  the  uniform  flow  is  not  predicted, 
while  in  the  downstream  regions,  additional  wake  effects 
have  not  been  taken  into  account.   However,  in  view  of  the 
large  area  of  good  agreement  with  experiment,  Wooler's 
distributed  vorticity  model  was  selected  as  the  foundation 
for  further  investigation.   It  is  apparent  that  this  model 
gives  a  fair  representation  of  the  actual  contra-rotating 
vortices  associated  with  the  jet-crossf low  interaction. 
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IV.  EXTENSION  OF  WOOLER  MODEL 

Upon  successful  confirmation  and  correction  of  Wooler's 
distributed  vortex  model  and  as  a  step  towards  developing  an 
improved  model,  it  was  felt  necessary  to  verify  the  trailing 
vortex  assumption,  that  these  vortices  are  tangential  to  the 
jet  rather  than  contained  in  the  jet.   Retaining  the  basic 
theoretical  concepts  for  calculating  the  incremental  circula- 
tion of  a  finite  length  of  the  jet,  containment  of  the  trailing 
vortices  was  approximated  by  successive  conforming  discretized 
vortex  elements.   These  vortex  elements  were  formed  by  restrict- 
ing the  length  of  the  trailing  vortices  to  the  linear  distance 
between  z,  and  z2,  the  arbitrarily-chosen  end  points  of  the 
finite  length  of  jet,  ds.   In  order  to  satisfy  the  Helmholtz 
vortex  theorem,  that  a  vortex  cannot  end  in  the  fluid,  the 
above  horseshoe  vortex  system  was  closed  by  connecting  the 
trailing  arms  with  another  bound,  contra-rotating  vortex. 
(See  Fig.  7)    By  making  this  alteration,  vortex  loops  of  in- 
creasing incremental  strength,  as  given  by  Eq.  (8) ,  can  be 
made  to  conform  to  the  curvature  of  the  jet,  as  in  the  actual 
flow  situation.   This  resulting  vortex  pattern  is  analogous 
to  that  created  by  a  continuously  accelerating  wing  /~~23_7. 
As  a  wing  accelerates,  the  strength  of  the  produced  vortices 
correspondingly  increase,  a  situation  very  similar  to  the 
change  in  elemental  vortex  strength  associated  with  the  jet 
curvature. 
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Using  this  contained-vortex  system,  the  interference 
velocity  due  to  the  jet  at  any  arbitrary  point  on  the  plate 
was  determined  by  applying  the  Biot-Savart  Law  (See 
Appendix  A).   In  this  application,  it  was  necessary  to 
approximate  the  trailing  vortices  as  straight-line  segments. 
However,  as  the  sizes  of  the  jet  elements  become  small,  the 
linear  approximation  can  be  made  to  conform  to  the  curva- 
ture with  an  exactness  that  is  only  limited  by  numerical  prac- 
ticalities.  The  total  interference  components  due  to  the 
entire  jet  result  from  the  summation  of  elemental  contribu- 
tions and  are  given  by 

2  N-> 


U    3  \U 


i=l 


i+1 


Ztan 

k=l 


.  /   Zk+1    zk 
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1  +  e 


k+1 


+  Z,J 


U    8  V  U  j 
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1  +  e 

where  N  is  the  number  of  jet  elements  and  z,  and  z,  , ,  are 

J  k      k+1 

the  endpoints  of  the  general  element.   Again,  the  parameters 
u,  v»  and  w  are  defined  as  geometric  coefficients  resulting 
from  applying  the  Biot-Savart  Law  around  the  element.   The 
subscripts  1,2,3  and  4  refer  to  the  forward  and  after  bound 
and  left  and  right  trailing  vortex  contributions,  respectively. 
These  are  given  by 
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where  the  various  £ '  s ,  ?'s  and  n's  are  given  by  Equations 
(Al,  A2,  A3)  of  ADpendix  A.   With  these  interference  velocity 
components  at  numerous  points  on  the  plate,  the  pressure 
coefficient,  Eq.  (9) ,  was  determined  for  selected  points 
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on  the  plate.   The  resulting  contours  of  constant  pressure 
coefficient  are  shown  in  Figs.  8,  9,  and  10. 

These  resulting  plots  show  a  somewhat  different  agree- 
ment with  the  experimental  data  /~~6_7  tnan  did  the  results 
of  Wooler's  model;  however,  the  discrepancy  can  be  directly 
attributed  to  the  manner  in  which  the  trailing  vortices  are 
treated.   In  Wooler's  model,  the  total  interference  velocity 
due  to  the  trailing  vortices  results  from  the  cumulative  el- 
emental vortices  acting  at  different  distances  relative  to 
a  point  on  the  surface,  while  in  the  conforming  vortex  model 
the  trailing  vortex  contribution  results  from  the  cumulative 
elemental  vortex  strength  acting  at  a  single  radius,  (See 
Fig.  11).   Therefore,  Cooler's  assumption  of  treating  the 

trailing  vortices  as  being  tangent  to  the  trajectory  intro- 
duced an  error  that  under  some  conditions  tended  to  improve 

agreement,  where  the  conforming  vortex  model  removed  this 
error  at  the  expense  of  some  loss  in  agreement,  at  least  at 
the  larger  values  of  U./tl.   However,  there  have  been 
studies  /~~12,2l7  that  present  the  concept  of  an  effective 
origin  of  the  jet  vortex  system  -  a  region  where  the  turbu- 
lent mixing  processes  extend  across  the  entire  jet  causing 
definite  deflection  of  the  jet  and  establishment  of  the  contra- 
rotating  vortices.   Schmidt  /~2l7  formulated  an  empirical 
expression  for  the  position  of  the  effective  source  as  a 
function  of  the  jet-to-mainstream  velocity  ratio  given  by 

Z  =  (15  a/o)"1 

where  a  =  U/U  .  .   Applying  this  effective  source  expression 
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to  the  conforming  vortex  model  produces  better  agreement 
with  experiment  (Fig.  12) . 

Comparison  of  the  conforming  vortex  model,  including  an 
effective  source,  with  experimental  data  shows  the  same, 
if  not  better,  far-field  agreement  with  an  improvement  in 
the  predicted  pressure  distribution  in  upstream  and  downstream 
regions.   Therefore,  this  model  which  more  closely  represents 
the  actual  flow  situation,  replaced  Wooler's  model  as  the 
basis  for  further  investigation  and  improvement. 
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V.  IMPROVEMENT  OF  CONFORMING  VORTEX  MODEL 

Evaluation  of  the  conforming  vortex  model  results  in- 
dicated the  continued  lack  of  agreement  with  experiment  in 
two  distinct  areas:   upstream  and  wake.   In  order  to  improve 
the  agreement  in  these  areas,  it  was  felt  that  additional 
corrections  could  be  formulated  from  a  qualitative  analysis 
of  the  flow  phenomenon.   The  discrepancies  in  the  upstream 
area  can  be,  at  least  partially,  attributed  to  the  blockage 
of  the  mainstream  flow  by  the  presence  of  the  jet,  while  the 
discrepancies  in  the  wake  area  are  caused  by  turbulent  flow 
separation,  vortex  shedding  and  turbulent  entrainment.   In 
view  of  the  present  uncertainty  concerning  the  exact  character 
of  the  wake  area,  it  was  decided  that  the  development  of 
corrections  for  the  wake  must  await  future  investigation  and 
experimentation.   Therefore,  the  upstream  blockage  correction 
was  undertaken. 

Near  the  point  of  injection,  the  physical  appearance  of 
a  circular  jet  injected  into  a  crossing  flow  is  that  of  a 
cylinder.   Interpreting  the  jet  as  a  cylinder  in  an  irrotational 
flow,  potential  flow  theory  predicts  perturbations  to  the 
uniform  flow  that  are  given  by 
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where  r  and  9  are  cylindrical  coordinates  of  points  in  the 
flow  field  and  a  is  the  radius  of  the  cylinder.   These 
velocity  perturbations  caused  by  blockage  can  be  combined 
with  the  interference  velocity  components  derived  from  the 
vorticity  formulation  to  adjust  the  coefficient  of  pressure 
(Eq.  9)  at  arbitrary  points  on  the  surface.   The  initial 
(and  most  simple)  assumption  was  to  make  the  diameter  of  the 
blockage  cylinder  the  same  as  the  diameter  of  the  jet  orifice 
The  resulting  pressure  distribution  indicated  an  excessive 
blockage  that  is  not  present  in  the  real  flow.   The  actual 
flow  does  not  produce  a  solid  interface,  as  would  a  cylinder, 
but  instead  an  entraining  interface  surrounding  a  core  of 
uniform  velocity  approximately  equal  to  U . ,  the  jet  exit 
velocity.   To  better  describe  this,  a  potential  cylinder  of 
diameter  less  than  the  diameter  of  the  jet  was  introduced. 
Using  many  combinations  of  cylinder  diameters  with  different 
jet-to-mainstream  velocity  ratios,  an  expression  for  a  in 
Eq.  (10)  was  formulated  in  terms  of  the  velocity  ratio. 

a  =  0.96  /a 

Applying  this  formulation  to  the  blockage-cylinder  perturba- 
tion velocity  calculations  results  in  pressure  contour  plots 
(Fig.  13,14,15)  that  show  improved  agreement  with  experiment 
in  the  far  field  and  upstream  areas  for  jet-to-mainstream 
velocity  ratios,  4  <  —  <  8.   However,  as  the  velocity  ratio 
increases,  the  model  results  and  experiment  diverge.   The 
apparent  cause  of  this  disagreement  can  be  resolved  by 
qualitatively  examining  the  actual  flow.   With  increasing 
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velocity  ratios,  the  entrainment  of  crossflow  fluid  by  the 
jet  becomes  the  more  dominant  interference  factor,  while  the 
blockage  factor  is  relatively  less  important.   Wu  and  Wright 
/T9  7  have  concluded  that  "the  amount  of  crosswind  fluid  en- 
trained increases  with  increasing  jet  speed"  when  the  cross- 
wind  speed  is  held  constant.   In  addition,  this  increased 
entrainment  appears  to  be  especially  important  in  the  down- 
steam  wake  region  and  causes  a  further  reduction  in  the 
blockage  effect  as  compared  to  that  due  to  a  solid  cylinder 
in  a  crossflow.   The  resulting  effect  of  the  increased  rate 
of  entrainment  at  higher  jet  speeds  is  larger  interference 
velocity  perturbations  in  the  flow  field  surrounding  the  jet 
orifice  near  the  plate.   In  view  of  this,  the  lack  of  agree- 
ment at  higher  velocity  ratios  can  be  explained  while  also 
Dointing  out  the  need  for  additional  correction (s)  to  account 
for  the  change  in  the  relative  importance  of  blockage  and 
entrainment  with  changes  in  velocity  ratio. 
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VI.  CURRENT  MODEL  CAPABILITY 

Realizing  the  limitations  of  the  current  analytical 
model,  it  is  nevertheless  appropriate  to  investigate  its 
predictive  capability.   Figure  16  shows  the  results  of  the 
analysis  expressed  as  the  force  induced  on  a  flat  plate  due 
to  a  jet  injected  at  right  angles  to  a  crossflow.   The  force 
is  expressed  as  a  fraction  of  the  ideal  jet  thrust  and  the 
negative  values  indicate  suction  forces  (in  a  direction  op- 
posite to  the  ideal  jet  thrust) .   Correlation  with  experiment 
/"~4_7  shows  a  difference  in  the  value  of  k  in  Eq.  (2)  .   The 
value  of  k  -  -4.4  resulting  from  this  analysis,  although  not 
in  the  predicted  range,  -3.0  <  k  <  -1.5,  still  displays  the 
correct  trend  and  provides  encouragement  as  to  the  potential 
success  of  a  model  based  upon  the  present  methods.   The  main 
difference  in  the  approximation  constant  results  from  the 
lack  of  proper  accounting  for  the  entire  entrainment  effect, 
a  factor  which  progressively  degrades  the  model  at  higher 
velocity  ratios.   However,  as  an  initial  design  approximation 
for  low  velocity  bow  thrusters  active  at  relatively  high 
ship  speeds,  the  conforming  vortex  model  presented  here  pro- 
vides a  satisfactory  prediction  of  the  interference  effects 
--  a  step  toward  improving  the  accuracy  of  lowering  and  con- 
trol determinations. 
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VII.  CONCLUSIONS 

The  conforming  vortex  model,  modified  with  corrections 
based  on  experimental  observations,  is  presented  as  an 
initial  formulation  of  a  simple  yet  physically  consistent 
representation  of  jet-crossf low  interference  effects.   Pres- 
sure distributions  in  the  far-field  are  adequately  predicted 
and  the  calculated  results  obtained  with  the  model  provide 
reasonably  good  agreement  with  experiment  for  the  lower 
jet-to-crossflow  velocity  ratios.   Agreement  is  less  than 
satisfactory  for  the  higher  velocity  ratios.   However,  the 
model  does  indicate  the  relative  importance  of  different 
flow  phenomena  with  changing  velocity  ratios.   The  most  evi- 
dent weakness  of  the  model  is  the  lack  of  an  adequate  method 
for  including  entrainment  and  blockage  effects  near  to  the 
origin  of  the  jet  at  higher  velocity  ratios.   By  restricting 
the  range  of  application,  the  conforming  vortex  model  demon- 
strates sufficient  agreement  to  be  utilized  as  a  first 
approximation  technique  for  designing  bow  thrusters. 
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VIII.  RECOMMENDATIONS 

Further  experiments  should  be  conducted  to  provide  more 
correlation  data,  especially  in  the  case  of  water  jet  into 
water  crossflow.   Empirical  constants  utilized  in  the  current 
model  and  those  that  will  become  necessary  when  trying  to 
model  increased  entrainment  can  be  refined  when  a  larger  data 
base  is  available.   The  encouraging  results  of  this  method 
of  analyzing  the  flow  field  should  not  be  overlooked,  and 
additional  corrections,  as  indicated  in  this  analysis,  should 
be  made  to  improve  its  accuracy. 
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Figure  1.   Sketch  of  bow  thruster-crossf low  interaction, 
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Figure  3.   Wooler  Horseshoe  Vortex  System 
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Figure  4.   Comparison  Of  Pressure  Coefficient  Contours 

Experiment  /  6_7>  Wooler  

Velocity  Ratio  =  4 
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Figure  5.   Comparison  Of  Pressure  Coefficient  Contours 
Experiment  /  6_7#  Wooler  


Uj/u  =  8 
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Fiqure  6. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment /  6  ~J ,   Wooler  

U./U  =  11.3 
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Figure  7.   Conforming  Vortex  Loop  System 
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Figure  8.   Comparison  Of  Pressure  Coefficient  Contours 

Experiment /  6  7 ' ,    Conforming  Vortex  Model 

U./U  =  4.0. 
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Figure  9.   Comparison  Of  Pressure  Coefficient  Contours 

Experiment  /  6_7,  Conforming  Vortex  Model 

U./U  =  8.0 
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Figure  10. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment  /  6_7*  Conforming  Vortex  Model 

U./U  =  11.3 
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Figure  a.   Wooler  Model 


Figure  b.   Conforming  Vortex  Model 


Figure  11. 


Comparison  Of  Wooler  (a)  And  Conforming  Vortex 
Model  (b)  Effective  Radii  For  Induced  Velocities 
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Figure  12. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment  /  6_7/  Conforming  Vortex  Model 

With  Effective  Source  Correction .U./U  =8.0 
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Figure  13. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment  /  6  7,    Fully  Corrected  Conforming 

Vortex  Model .U./u  =4.0 
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Figure  14. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment /  6  7 ' ,    Fully  Corrected  Conforming 

Vortex  Model — .  U.7u  =8.0 
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Figure  15. 


Comparison  Of  Pressure  Coefficient  Contours 

Experiment /  6  7 ,    Fully  Corrected  Conforming 

Vortex  Model ,0  .7v   =  11.3 
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Figure  16. 


Variation  In  The  Induced  Force  Due  To  Jet-Crossf low 
Interaction 
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APPENDIX  A 
FORMULATION  OF  INTERFERENCE  VELOCITY  COMPONENTS 

The  formulation  of  the  interference  velocity  due  to  the 
vortex  loop  system  of  the  conforming  vortex  model  is  given 
below.   The  formulation  requires  the  use  of  three  separate 
coordinate  systems  that  are  functions  of  the  endpoints  of  a 
jet  element  and  the  position  of  an  arbitrary  point  on  the 
plate.   From  Fig.  (Al) ,  these  are  given  by: 


5-j   «  -  |z,sina,  +  (x,  -  x  )  cosa,l 


(Al) 


n  =  y 


5   =   z,cosa,  -  (x,  -  x  )  sina, 
*1      1    11    p      1 

for  the  upstream  bound  vortex 

C2   =     z2sina0  +  (x2  -  x  )  cosa2|  (A2) 

n  =  y 
C~   =   z_,cosa-  -  (x-  -  x  )  sina. 

2.  2  2  2  p        2 

for   the   downstream  bound  vortex,    and 
*Bi    =  "  [ZisinaB    +     (Xi    "    V     COSaB] 

n  =  yp  (A3) 

CB  =   z.cosaB  -  (x.  -  xo)  sinaB 
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for  the  trailing  vortices,  where  x   and  y   are  the 

P      P 

coordinates  of  an  arbitrary  point  on  the  plate  and  x, ,  z, 
x2,z2  are  the  coordinates  of  the  endpoints  of  an  element 
of  jet.   The  angles  a, ,  a~  and  a_  are  the  angles  between  the 
x-direction  and  respective   £'s  and  are  given  by 

-1 
a,      =      tan 


[csch    (z,/B)l 
|csch    (z2/B)| 
aB     =      tan~         (z2   -    z1)/(x2   -   x^J 


.   -1 
a,   =   tan 


To  determine  the  induced  velocity  at  an  arbitrary  point, 
the  law  of  Biot  and  Savart 

L 
must  be  individually  applied  to  both  the  bound  and  trailing 
vorticies.     Due  to  similarities  in  the  integrations,  only 
one  bound  and  one  trailing  vortex  formulation  will  be  pre- 
sented and  necessary  changes  applied  for  the  other  vortices. 
The  circulation  of  a  finite  amount  of  jet,  given  by  /~~Eq. 
(  8)7 


k    =   J-  ua(!i)2 


,   z2/Bd        ,   z./Bd 
tan   e      -  tan  J  e 


is  the  incremental  contribution  which  is  added  to  each 
successive  element,  such  that  the  total  circulation  for 
a  given  element  is  given  by 
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rn     '     }™W 


n 


,    -i    zk+i/Bd     ,    -i    VBd 

tan        e  -   tan        e 


(A5) 


Therefore,  the  induced  velocity  due  to  any  of  the  vortex 
filaments  associated  with  the  n-th  element,  1  £  n  £  N, 
where  N  is  the  number  of  jet  elements,  can  be  determined 
from 

t   i     n  /  sin$  ,.,.. 

Nnl   =  4?  J  ""I"  ds  (A6) 

-'   r 

L 

where  L  is  the  length  of  the  filament,  3  is  the  angle 
between  the  vortex  filament  and  the  radius,  r,  to  an 
arbitrary  point  on  the  plate. 

For  the  forward  bound  vortex  (See  Fig.  (A2)) 

sm3  =  s 


v«i2  +  h 


2  2 

+  (n  +  s)4 


and 


2        ^2       2  x2 

r     =   5X   +  Z,1     +    (n  +  s) 

Substituting  these  equations  into  Eq.  (A6)  gives 

rn  r        V*!2  ♦  cx2' 

lqnl'  =  4?     T — 2 2 rf  3/2  ds       (A7) 

L 
the  magnitude  of  the  induced  velocity  at  an  arbitrary  point 

on  the  plate  due  to  the  bound  vortex  of  the  n-th  element. 
The  x-direction  component  of  this  velocity  is  given  by 
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qxl     =    'qnl'      cosy^osc^  -   siny-j^sinaJ 


(A8) 


where   cosy,    = 


1  VI7^7 


and   sinYi    = 


vim? 


Inserting  Eq.  (A7)  into  Eq.  (A8)  and  applying  the  integration 
limits  gives 

r  Ji 

V 

_        n 

qxl     Ttt 


Cx2  +  *>*   +  (n  +  s)2 


_  v jjy   ds  cosa,  - 


(A9) 


-k 


2      2  2 

Cl   +  ?1   +  (n  +  S) 


1372"  dS  Sinal 


Integration  of  Eq.  (A9)  results  in 


Lxl 


n 

47T 


n  +  h 


V  ^i   kL2  +  t*  - 


-  h 


\ 


cosa. 


n  +  h)2  yk\   +  £2  +(r\-H)7 


X2  r2 


n  +  % 


+  Vs 


i2  +  C2  +  (n  +  ^s)2 


S2  +  C2  +  (n-^)2 


■sina. 


(A10) 

the  x-direction  component  of  induced  velocity  due  to  the 
forward  bound  vortex.   Similarly,  the  x-direction  component 
of  induced  velocity  due  to  the  downstream  bound  vortex  is 
given  by 
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lx2   "   Ttt 


n  +  H 


n  -  % 


cosa. 


C2     +  ^2   Vv422  +  ?22  +    (n  +  ^)2    Vc22+c22+(n-%)27 


n  +  ^ 


n  -  % 


C2    +  c2  ^22+^22+(n^)2'   V?22+c22+(n^)2y 


sina. 


(All) 


The  approximation  of  the  bound  vortices  as  being  linear 
eliminates  the  y-direction  component  of  induced  velocity 
and  since  an  image  system  is  utilized  to  formulate  the  flat 
plate,  it  is  unnecessary  to  calculate  any  z-direction 
component  of  induced  velocity.  - 

For  the  left  trailing  vortex  filament,  as  viewed  from 
upstream  (See  Fig.  (A3) ) . 


sing  = 


v^T 


(n  -  H) 


V^7 


2     2 

(  n-  h)     +  sz 


2    r   2    /     M2     2 
r=CB+(n-is)   +s 


Insertinq  these  expressions  into  Eq.  (A6)  gives 

r 


ln3 


n 
4iT 


SB2  +  (n  -  h)2 


[cB2  +  (  n  -  H)2   +  s2] 


37? 


ds 


(A12) 


the  magnitude  of  the  induced  velocity  at  an  arbitrary 

point  on  the  plate  due  to  the  left  trailing  vortex  of  the 

n-th  element.   Resolving  this  velocity  into  x  and  y  com- 
ponents gives 
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Lx3 


ln3 


cosy.  sinaD 


and 


(A13) 


lY3 


ln3 


smy 


B 


where  cosyB  = 


n  -  h 


?B2  +  (n  -  h)2 


and  siny_  = 


■B 


CB2  +  (n  -  h)2 


Inserting  Eq.  (A12)  into  Eq.  (A13)  and  applying  integration 
limits  gives 


r      r  B2  i. 

!x3   "         4tt     /        \7T~.      7 ,2^    213/2 
h^B   +<*-*>    +S    J 


ds    sina 


B 


r        r    B2  c„ 

-     n      I  B  ^1 

qy3      1*     /     rT2T777r72772T372"  c 


F    Ic  2+(n-is)2+s 

Integration  of  Eq.  (A14)  yields 


(A14) 


»x3 


n 
4it 


n  ■  H 

L^B2  +  (n  -  h) 


'32 


•Bl 


/cB2+(n^)2+C22V^B2+(n^)2+C2{ 


sina. 


n 
qy3  *  47 


'B 


'B2 


'Bl 


LV+^^^W^"2  VcB2+(n-*)2+^l' 


(A15) 

the  x  and  y  comoonents  of  induced  velocity  due  to  the  left 
trailing  vortex.   Similarly  the  x  and  y  components  of  the 
induced  velocity  due  to  the  right  trailing  vortex  are  given 
by 
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In   D  +  * 

qx4  '  4tt   ~T~ 


'B2 


•Bl 


<b    +  ">  +  *>    \VcBa+(n*wa*«2a   V^W^r 


sma. 


n 


ly4    4tt 


'B2 


'Bl 


L^  +  <i  +  *>2  V,B2+(n^)2+4  V.B2+  (n^)2+^Jj 


(A16) 
Therefore,  the  total  x  and  y  components  of  induced  velocity 
due  to  an  entire  element  of  jet  are 


qxT  =  qxl  +  qx2  +  qx3  +  qx4 


qyT  =  qy3  +  qy4 


(A17) 


Inserting  Eqs. (A10) ,  (All),  (A15) ,  and  (A16)  into  Eq.  (A14) 
yields 

rn 

qxT   =  Tn    LU1   cosctl    "    u2    cosa2~wisinai+w2sina2+{"w3+w4)sinaBj 


n 
qyT        Ttt 

where 

u,    = 


V3    "    V4 


3 


(A18) 


n  +  h 


n  -  h 


Ci    +  ci    LVc12+C12+(n+!s)2    Vc12+?12+(n-^)2. 


U2   =   —    .     .2 


n  +  h 


n  -  Jj 


?2     +  ?2   LVs22+c22+(n+%)2    V522+c22+(n-%>2- 


V3   =  71 


'B 


v.  »  — 7 


cB    +  (n-*t) 


'B 


'B2 


'Bl 


-VejL+c2  +  (n^)2    Vd,+?2  +(n-%)2- 


*B2    ^B 


Bl    SB 


'B2 


'Bl 


?b    +   <"♦■.!    LVd,+CB2+(n^)2    Vd,+?02+<n^)2J 


'B2    *B 


Bl    ^B 
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wl  H  ,  2  .  r  2 


n  +  J5 


n  -  *s 


w2    =   — 1 2 

^2  g2 


'-Vc12+;12+(n+*3)2     V?12+C12+(Ti-i5) 
n  +  h  r\  -  H 


•V?22+c22+(n+^)2    V?22+c22+(n-is)2-' 


W3    =    — Z 


n  -  % 


cB     +    (n-*j) 


•B2 


*B1 


L^b2+^b2+(^)2     V5^+cB2+(n-i|) 


W4  =  71 


n  +  H 


CB     +    (n+Js) 


'B2 


•bi 


^B2+?B2+(n+!s)2  V5Bi+5B2+(n+h)2j 


Substituting  the  expression  for  T   from  Eq. (A5)  into  Eqs. 
(A18)  and  non-dimensionalizing  yields 


U    8  V  U 


2    N 


i=H 


-w,  sin  a.  +w-sinct.  , ,  +  (-w,+w. )  sin  a _+u,  cos  a.  • 
1    l   2     l+l    3   4     B   1    l 


u2cosa.+1 


x  )   tan 


zk+l    zk 
s      -  e 


k=l 


1  +  e 


(zk+l+V 


N 


5-UW     Eh-M 


tan 


i=l 


k=l 


zk+l    zk 
-1  /  e -  e 

"  7  (zk+i+zk)- 

1  +  e 
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Figure  Al.   Coordinate  Systems  For  Conforming  Vortex  Model 
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Figure  A2.   Geometry  Of  Bound  Vortex  Integration 
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*  p 


n-*5 


-  ~-^ 


Figure  A3.   Geometry  Of  Trailing  Vortex  Integration 
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xxxxx>xx>>xxxxxxxxxxxxxxxxxx 


X 
X 
X 
X 
X 
X 
X 
X 
X 


THIS  PROGRAM  SIMULATES  A 
FLUID  JET  I  MECTEC  INTO 
A  LMFGPM  CROSSFLCW  AND 
FURTHER  DETERMINES  THE 
PRESSURE  OISTRI8LTICN  QN 
THE  SURFACE  PROM  1*HICH 
THE  JET  IS  EJECTEC 


XXXXX>XXXXXXXXXXXXXXX>XXXXXX 


DEFINITION  CF  TERMS 


C  NTYFE   -  MODEL  TYPE  :  C-CCNFORM 

C  1-CCNFORM  8LGCKAGE 

C  2-RAISED  CONFORM 

C  3-RAISED  CONFCRN  8LCCKAGE 

C  (INPUT) 

C 

C  NRAC    -  NUMBER  OP  RACIAL  PGSITICNS  CN  SURFACE 
C  (INPUT) 

C  NANG    -  NUM8ER  OF  ANGULAR  POSITICNS  ON  SURFACE 
C  (INPUT) 

C  NPCINT  -  NUMBER  OF  ARBITRARY  POINTS  CN  SURFACE 
C  (NRAD*NANG) 

C  LJ        JET  VELOCITY  (INPUT) 

C  LV        MAIN  STREAM  VELOCITY  (INPUT) 

C  C         ACTUAL  DIAMETER  OF  JET  CRIFICE  (INPLT) 

C  P         ISON-CIMENSICNAL  RAOIAL  PCSITION  OF  SURFACE 
C  PCINTS  (DIMENSIONED  AT  LEAST  NRAD)  (INPUT) 

C  ThETA   -  ANGULAR  PGSITICN  OF  SURFACE  POINTS 
C  (DIMENSIONED  AT  LEAST  NANG) 

C  >P        X-POSITION  CF  SURFACE  POINTS 
C  (CI^ENSIONED  AT  LEAST  NPCINT) 

C  VP        Y-POSITICN  CF  SURFACE  PCIMS 
C  (DIMENSICNEC  AT  LEAST  NPCINT) 

C  X         HCRI20MTAL  COCRDINATE  OF  JET 

C  2         VERTICAL  CCCRCINATE  OF  JET 

C  2INC    -  NON-CIMENSICNAL  INCREMENTAL  STEP  SIZE  FCR  Z 

C  CELTAZ  -  2-POSITION  CF  EFFECTIVE/VIRTUAL  SOURCE 

C  CELCCK  -  EFFECTIVE  BLCCKAGE  CYLINCEP  CIAMETEP 

C  F.KSAV   -  CUMULATIVE  VCRTICITY  ALCNG  THE  JET 

C  XI        NATURAL  CCCRCINATE  TANGENT  TO  JET 

C  iET*    -  NATURAL  CCCRCINATE  NORMAL  TC  JET 

C  ETA     -  NATURAL  COCRDINATE  P ERPENC  ICULAR  TC  >I 
C  AND  2ETA 

C  ALPHA   -  ANGLE  BETWEEN  X-AXIS  AND  XI-AXIS 

C  FAREA      -    INCREMENTAL    PLATE    AREA 
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C      TCTFCR  -  TOTAL  RESULTANT  FORCE  ON  PLATE 
C      >eAG    -  X-PCSITION  CP  RESULTANT  FCFCE 
C      YT        NOPMJLIZEC  INCUCEO  FORCE 
IMPLICIT  R£AL*8  <A-h,C-Z) 

CINENSICN  U(5C,50)  ,V ( 50 ,50 ) , M <50 , 5C J ,R ( 50 ) , TF ETA ( 50  )  t 
X  CP(50,50) 

CI  PENSION  XP<450) ,YF(45C),ZP(450) 
C 

C   CAT*  INPUT 
C 

PEAC(2,60)  NRAD 

PEAC<2,60)  NANG 

C    ZEPC  CLT  STORAGE  MATRICES 
f^PCINT  =  NRAO*NANG 

CC  1C  I  =  ifNFOINT 

XP(  I)  =  0.0 
VPU)  =0.0 

ZP  <  I  )   =  0.0 

1C     CONTINUE 

CC  30  I  =  1,NPAD 

R(I)  =  0.0 

CO  2C  J  =  1,NANG 

THETA(J)  =  O.C 

UI  tJ)  =  0.0 

V(  I, J)  =  CC 

CP(ItJ)  =  0.0 
20     CONTINUE 
30     CGNTINUE 

REAC(2t60}  NTYPE 

REAC(2,40)  UJ,UN,D 
k»RITE(6,50J  U^,UM,0 

REAC(2,70J  ZINC 

REAC(2,70)(R<  I)  ,I=1,NRAC) 


40 
5C 

6C 
<C 

8C 

SC 

c 
c 

C 


VELOCITY    RATIO 

VR    *    LJ/UM 
FCSITICN    OF    EFFECTIVE    SCUPCE 

CELTAZ    =    VR*DSQRT(VP)/15.0 
CCRRECTICN    FCR    DELTAZ    =    0.0 

IFUTYPE.LE.l)    CELTA2    =    0.001 

kRITE(6t8C)    OELTAZ 
fcRITE<6,90)     ZINC 

FCRNAT(3F8.4) 

F0RNAT(,0,»2X,f JET    VELOCITY    =    ',     F8  .4 ,2X , • UN  I  FCRN ■  , 
X'VELCCITY    =    •  ,F8.4t2X,  'JET    CIAMETEP    '    »,F8.4) 

FCRMTU4) 

FCR*AT(F1C5) 

FCPM4T(,0,t2Xi 'EFFECTIVE  SOURCE  IS,iF5.2,2XT 
>*CIANETERS  AeCVE  THE  PLATE') 

FORVAT  (,C,f2>,'Z  STEP  SIZE  =  »tF5.3) 

CEFIMTION  CF  PI 
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PI    =    2.0*ARSIM1.00CCCQ0O) 

TRJJECTCRY    CONSTANT 

E    =    C.1<3*VR**2 

CALCULATION  OF  INTERFERENCE  VELOCITY 

L  =  C 

CALCON  =  DFLOAT(NANG)  -  1.0 

C    INCREMENT  RADIUS  (INPLT  VALUES) 

CC  140  I  -    1,I*RAD 
RCALC  =  R(I ) 

C    INCREMENT  ThETA  (5  OEGREE  INCREMENTS) 

CO  130  J  =  1VMNG 

K  =  J-l 

THETA(J)  -    DFLQAT(K)*  (PI/CALCON) 

C    SMFT  TO  CARTESIAN  CCCPCINATES 

L    =    L-H 

COST    =    DCOS(TFETA( J)) 
SINT    =    OSIN(THETA( J)) 
XP(L)    =    RCALC*COST 
VP(L)    =    RCALC*SINT 
C 

21    =    0,0 

11    =    C.00CC1 

HKSAV    =    0.0 

ICC         21    =    11 

11    *         Zl    ♦    ZINC 
211    -    Zl    ♦    OELTAZ 
Z22    =    Z2    +    DELTAZ 

>1    =    B*(DC0SH<21/3)-1.0) 
>2    =    B*(DCOSH(Z2/B)-1.0) 

TANAE  =  ( Z2-Z1)/(X2-X1) 
TANA  1  =  1.0/DSINHC Zl/B) 
TANA2  =  1  .O/DSINHC  22/e) 

/SLFCNE  =  CAT4MTANA1  ) 
4LFTW0  =  CATAMTANA2) 
ALFEAR  =  DATAN(TANAE) 

HNM  '  DSIMALFONE) 
SINA2  =  DSIN(/LFTW0) 
SIN/E  =  CSIN<4LFBAR) 

C0S41  =  CCCS(ALFONE) 
CCSA2  =  0C0S<ALFTW0) 
C0SA6  =  OCGS(ALFBAR) 

C   NATLRAL  COORDINATES  IN  TERPS  OF  CARTESIAN  COORDINATES 

XICNE  =  -(Z1HSINA1  +  (Xl-XP(L) )*CCSA1> 
XITVC  =  -(Z22+SINA2  «•  ( X2-XP< L ) ) *CCS A2 ) 
XIB4R1  =  -(Z11*SINAE  +  (Xl-XP(L)  MCCSAB) 
XIE4R2  =  -<222*SINAE  +  < X2-XP (L ) ) *CCSAB) 

ETA  =  YP(L) 

2ET41  =  (Z11*C0SA1  -  (X1-XP<L))*SINA1) 
ZETA2  =  <Z22*C0SA2  -  ( X2-XP ( L ) ) *S IN IZ ) 
2ETA6    =    <2U*CCSAB    -    (  Xl-XP  (L  ) )  *S  INAE  ) 
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C   GECNETPIC  CONSTAMS 

ETPL  =  ETA  4  C.5 

ETNI  =  ETA  -  C.5 

SA  =  XICNE**2 

S6  =  XITW0**2 

SC    =  XIBAR1**2 

SC  =  XIBAR2**2 

SS  -    ZETA1**2 

$F  =  ZETA2**2 

SG  -  ZETAE**2 

SF  =  ETPL**2 

SI  '    ETMI**2 

CA  =  DSCRTiSA  4  SE  4  SH) 

C8  =  OSQRTtSA  +  SE  4  SI  ) 

CC  =  DSCRKSe  4  SF  4  SH) 

CE    =    CSQRT(S8    4    SF    4    SI  ) 

CF    =    CSQRKSD    +    SG    •»    SI  ) 

£G    =    CSGRKSC    4    SG    +    SI) 

CH    =    DSQRT(SC    +    SG    +    SH ) 

CI    =    CSGRKSC    t    SG    ♦    SH) 

CCM    =    <ETPL/CA    -    5TMI/CB)/(SA    4    SE) 
CCN2    =    (ETPL/CC   -    ETMI/CE)/(S8    4    $F) 
CCIS2    =    (XIBAR2/OF    -    XI EAR1/CG)/ < SG    <    SI) 
CGN4    =    (XIBAR2/DH    -    XIB AR1/DI ) / ( SG    4    SH) 

LI  =  Z5TA1*CCM 

L2  =  ZETA2*CCfv2 

Vi  *  ZETAe*CCN3 

V4  =  ZETAB*CON4 

V»l  =  XICNE*CGM 

V>2  -  XITWC*CON2 

V3  =  ETMI*CGN2 

fc4  =  ETPL*CCN4 

C       INCREMENTAL    VORTICITY    TERM 

FF    -     ((DEXP(Z2/B)-DEXP(Z1/B))/<1.C4CEXP((Z24Z1)/E))  ) 
FK    =    U. 0/8*0)*    VR**2    *    DATAN(FF) 
FK    =    HK    4    HKSAV 
FKSAV    =    HK 
INCREMENTAL     INTERFERENCE    VELOCITY    COMPONENTS 

LP    =    CO 

LP    =    FK*U-W1*SINA14W2*SINA24(-W34W4)*SINAB4L1*CCSA1  ) 

X    -<L2*COSA2M 


C 


\«P    =    0.0 

\P    =      HK*<V3    -    V4) 


L(I,J)    =    U(I. J)    4    UP 
WI  ,J)    =    V(  I, J)    4    VP 


TEST  FCP  CONVERGENCE 

TEST  =  CABS(UF)  4  DABS(VP) 

IF  (TEST.LE.C.COO0C1)  GC  TO  110 
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GC  1C  100 
C 
110    CONTINUE 

C   CALCULATION  OF  BLOCKAGE  VELOCITY  TERN 

IF(NTYPE.EQ.C)  GO  TC  120 

IF<MYPE.EC2  )    GO    TC    120 
C 

PeLCCK    =    C96/CSQRT(VR  ) 

CBLCCK    =    2.0*RBL0CK*C 

PK    =     (<R8LCCK/R< I)  )**2)/2.0 
LC    =    RK*(SINT**2    -    CCST**2) 
VC    =    -2.0*RK*SINT*CCST 
L(  I, J)    =    U<  It  J)    +    UC 
V(I  tJ)    =    V<  I* J)    +    VC 

C       CALCULATION    CF    PRESSURE    COEFFICIENT 

120  CPUtJ)     =   -4.C*(U(ItJ)    +    U(I,J)**2    *    VU,J)**2) 

C        CHANGE  CESIGNATICN  OF  VARIAELE 
C 

2P(L)  =  CFdt  J) 


C 

130 
140 


CONTINUE 
CONTINUE 


C   CALCLLATION  OF  INCREMENTAL  FORC? 

RG  =  0.0 
RI  =  0.0 
AREA  =  0.0 
PAREA  =  0.0 
FCPCE  =  0.0 
YNCI*  =  0.0 
TCTFCR  =  0.0 
TYNCN  =  0.0 

NAREA  -    NPAC  -  1 

CC  17C  I=1,NAFEA 
K  =  I  +     1 
PI  =  R(I> 

RG  -    R(K) 

C   CALCULATION  CF  IhCREMENTAL  AREA 


AREA  =  PI*(R0**2  -  PI**2) 
PAREA  =  AFEA/ <2.0*CALC0N) 
CO  160  J=1,NANG 

C   CALCULATION  CF  AVERAGE  PRESSURE  COEFFICIENT 

CPAVG  =  (CP(ItJ)  ♦  CP(K,jn/2.0 
C   CALCULATION  OF  EFFECTIVE  FCRCE 
FORCE  =  PAREA*CPAVG 

IF(v.EQ.l)  GC  TO  15C 

IF(J.EQ.NANG)  GO  .TC  150 
FCRCE  =  2.0*FCRCE 

C   CALCULATION  CF  MOMENT  ARM 
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15C    *P    =  UPC  +  RI)/2.C)*CCCS(THETA(  J)  ) 
C   CALCULATION  CF  MCMENT  ABOUT  Y-AXIS 

^NC>  =  FORCE*XM 
C   SUNNATICN  CF  T0T4L  FORCE 

7CTFCR  -  TOTFCR  +  FCRCE 
C   SUNNATIQN  CF  TOTAL  MOMENT 

TYPCV  =  TYMCN  +  YMCN 

160    CCNTINUE 
IK    CONTINUE 

C   CALCULATION  CF  CENTER  CF  ACTION 

>E AR  =  TYMCM/TOTFOR 
C   CALCULATION  CF  NCRMALIZEC  INDUCED  FORCE:  Y/T 

\RIIV\  =  1.0/VP 

YT  =  2.0*TCTFCR*VRINV**2/PI 

C   CAT/S  CUTPUT 

V»RI7fc<6,260)  C8L0CK   ' 
bRITE(6,220)  TCTF0R,XE4P 
V»RI1E(6f240)  YT,  VRINV 

180    CC  220  I  =  1,NANG 

ANGLE  =  ThETA(I>*18C.0/PI 

fcRI7E(6,180)  ANGLE 
19C    FCPVAT(«0« f2X  ,»THETA  =  '  ,F6  .2 t 8X ,  •  P • ,9X, • CP * ♦ «X » ■ L» 

x  9x,«v« ) 

CC  210  J  =  lthRAO 

VRITE(6,2CCJ  P(J),CF(J,I),U(J,I),VU,I) 
2CC    FCRVAT(»  ■ ,23X,F6.2,2X,F8.2,3X,F6.2,4X,F6.2) 
210    CGNTINUE 
22C    CONTINUE 

230         FCRMTCO'  ,2X,  •RESULTANT    FORCE    OF    •  ,  IX,  F7.2,  IX  , 
X'KTS    AT    X    =    «  ,F6.2  ) 

240         FCRMT(«0',2X,»Y/T    =    » tF6.2 ,1 X,  'FOR    UM/UJ    =     »,F6.3) 

^RITE(7,25C)(XP(I)»YP(I),ZP(I),I=l,hFCINT) 
250         F0RNAT(3F20.4) 
260         FCFNATC0S2X, 'EFFECTIVE    BUCCKAGE    OLINDER1, 

X'CIANETER    =    •  ,F4.2) 
27C  STOP 

ENC 
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